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We present a simple point process model of 1//^ noise, covering different values of the exponent 
f3. The signal of the model consists of pulses or events. The interpulse, interevent, interarrival, 
recurrence or waiting times of the signal are described by the general Langevin equation with the 
multiplicative noise and stochastically diffuse in some interval resulting in the power-law distribu- 
tion. Our model is free from the requirement of a wide distribution of relaxation times and from 
the power-law forms of the pulses. It contains only one relaxation rate and yields 1//^ spectra 
in a wide range of frequency. We obtain explicit expressions for the power spectra and present 
numerical illustrations of the model. Further we analyze the relation of the point process model of 
1/f noise with the Bernamont-Surdin-McWhorter model, representing the signals as a sum of the 
uncorrelated components. We show that the point process model is complementary to the model 
based on the sum of signals with a wide-range distribution of the relaxation times. In contrast to 
the Gaussian distribution of the signal intensity of the sum of the uncorrelated components, the 
point process exhibits asymptotically a power-law distribution of the signal intensity. The devel- 
oped multiplicative point process model of l/f@ noise may be used for modeling and analysis of 
stochastic processes in different systems with the power-law distribution of the intensity of pulsing 
signals. 

PACS numbers: 05.40.-a, 72.70.+m, 89.75.Da 
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I. INTRODUCTION 



1/f fluctuations are widely found in nature, i.e., the power spectra S(f) of a large variety of 
physical, biological, geophysical, traffic ^financial and other systems at low frequencies / have 
1/fP (with 0.5 < (3 < 1.5) behavior □ ODD . Widespread occurrence of signals exhibiting 
such a behavior suggests that a generic mathematical explanation of 1/f noise might exist. 
The generic origins of two popular noises: white noise (no correlation in time, S(f) ~ 1//°) 
and Brownian noise (no correlation between increments, S(f) ~ 1/f 2 ) are very well known 
and understood. It should be noted, that the Brownian motion is the integral of white noise 
and that operation of integration of the signal increases the exponent by 2 while the inverse 
operation of differentiation decreases it by 2. Therefore, 1/f noise can not be obtained 
by simple procedure of integration or differentiation of such convenient signals. Moreover, 
there are no simple, even linear stochastic differential equations generating signals with 1/f 
noise. Recently we derive a stochastic nonlinear differential equation for the signal exhibiting 
1/f noise in any desirably wide range of frequency j^. The physical interpretation of this 
highly nonlinear equation is not so clear and straightforward as that of the linear Langevin 
equation, generating the Brownian motion of the signal with 1/f 2 spectrum. Therefore, 1/f 
noise is often represented as a sum of independent Lorentzian spectra with a wide range 
of relaxation times Summation or integration of the Lorentzians with the appropriate 
weights may yield 1/f noise. 

Not long ago a simple analytically solvable model of 1/f noise has been proposed 0, 
analyzed & an d generalized [l(3]. The signal in the model consists of pulses or series 
of events (a point process). The interpulse times of the signal stochastically diffuse about 
some average value. The process may be described by an autoregressive iteration with a 
very small relaxation. The proposed model reveals one of the possible origins of 1/f noise, 
i.e., random increments of the time interval between the pulses (the Brownian motion in the 
time axis), sometimes resulting in the clustering of the signal pulses 

The power spectral density of such point process may be expressed as 

S(f)~2PfP k (0)/f. (1) 

Here f = (r^) is the expectation of the interpulse time r k = t k+ i — t k , with {t k } being the 
sequence of the pulses occurrence times or arrival times t k , whereas P k ij k ) is a steady state 
distribution density of the interpulse time r k in fc-space and I is the average intensity of the 
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signal 



I(t) = ^A k (t-t k ). 



(2) 



Function A k (t — t k ) represents the shape of the fc-pulse of the signal in the region of the 
pulse occurrence time tk- 

It is easy to show that the fluctuations and shapes of Ak (t — tk) for sharp pulses mainly 
influence the high frequency power spectral density. Therefore, in a low frequency region we 
can restrict our analysis to the noise originated from the correlations between the occurrence 
times tk- Then we can simplify the signal to the point process 



j (t) = a $ (t - t h 



(3) 



with a being an average contribution to the signal of one pulse or one particle when it crosses 
the section of observation. 

Point processes arise in different fields, such as physics, economics, cosmology, ecology, 
neurology, seismology, traffic flow, signaling and telecom networks, audio streams, and In- 
ternet (see, e.g., 0, 111 12, 111 3 an d references herein). The proposed point process model 
3, 0, Q] can been modified and useful for the modeling and analysis of self-organized sys- 
;ems jig , atmospheric variability , larg e flares from Gamma-ray Repeaters in astronomy 
~y\ . particles moving in viscous fluid jl8|, d ynamical percolation 1/f noise observed 

financial markets (21 1. cognitive experiments 



in cortical neurons and earthquake data 
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22j . the Parkinsonian tremors j23j, and time intervals production in tapping and oscilla- 
tory motion of the hand 2^ |. 

The analytically solvable model and its generalizations ll(j ] contain, however, 



some shortage of generality, i.e., it results only in exact 1/f (with (3 = 1) noise and only if 
P k (rk) — const when r k —>■ 0. On the other hand, the numerical analysis of the generalized 
model with different restrictions for diffusion of the interpulse time T k reveals 1/f^ spectra 
with 1<(3< 1.5 Q- 

The aims of this paper are to generalize the analytical model seeking to define the variety 
of time series exhibiting the power spectral density S(f) ~ 1/f^ with 0.5 < (3 < 2 and 
to analyze the relation of the point process model with the Bernamont-Surdin-McWhorter 
model p , representing the signal as a sum of the appropriate signals with the different rates 
of the linear relaxation. 
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II. POWER SPECTRAL DENSITY OF THE POINT PROCESS 



The point process is primarily and basically defined by the occurrence times t k . The power 
spectral density of the point process (j3J) may be expressed as 



S(f) 
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where T = tr — t«^> uj Ms the observation time, u = 2irf, and 



A (k; q) = t k+q - t k 



k+q-l 

£<■< 



(5) 



is the difference between the pulses occurrence times t k+q and t k . Here fc min and /c max are 
minimal and maximal values of index k in the interval of observation T and the brackets 
(...) denote the averaging over realizations of the process. 

It should be stressed that the spectrum is related to the underlying process and not to a 
realization of the process j3, [^J • Therefore, the averaging over realizations of the process 
is essential. Without the averaging over the realizations we obtain the squared modulus 
of the Fourier transform of the data, i.e., the periodogram which is fluctuating wildly and 
its variance is almost independent of T |25l. l26j|. For calculation of the power spectrum of 
the actual signal without the averaging over the realizations one should use the well-known 



procedures of the smoothing for spectral estimations 
Equation (J3J) may be rewritten as 
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/ . _2 N fcmax-g 

S (/) = 2a 2 u + lim ( -A £ £ cos [a; A (k; q)} 



where N = k max — fc min and 
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T 



(6) 



Q — 1 — &min 



lim 

T^oo 



N+l 

T 



is the mean number of pulses per unit time. The first term in the right-hand-site of Eq. (JBJ) 
represents the shot noise, 

Sshot (/) = 2a 2 u = 2dl, (7) 
with I — av being the average signal. 



Eqs. P))-(j7|) may be modified as 



N / I l\ 

S (f) = 2a 2 £ h ~ Y) Xa ^ ^ ® 

q=-N ^ ' 

and used for evaluation of the power spectral density of the non-stationary process or for 
the process of finite duration, as well. Here 

+00 



Xm H = ^ Hq) ) = J e lujA{q) ^ q (A (q)) dA (q) (9) 

—00 

is the characteristic function of the distribution density \l/ g (A(g)) of A(g), a definition 



A (g) = —A (— g) = A (k; q) is introduced, and the brackets (. . .) denote the averaging over 
realizations of the process and over the time (index k) P, . For the non-stationary process 
or process of the finite duration one should use the real distribution ^> q (A (q)) with the finite 
interval of the variation of A (q) or calculate the power spectra directly according to Eq. 
0. 

When the second sum of Eq. (jHl) in the limit T — > 00, due to the decrease of the 
characteristic function XA(g) (^) f° r finite uj and large q, approaches to zero, 

1 N 

lim - lei XA( q ) (w) -> 0, 

1 — >oo ± — 

q=~N 

we have from Eq. (jSJ) the power spectrum in the form 

= & (^E ei " A(fe;9) ) = 2j2f E (io) 

°° \ fc,g / g=-JV 

III. STOCHASTIC MULTIPLICATIVE POINT PROCESS 

According to the above analysis, the power spectrum of the point process signal is completely 
described by the set of the interpulse intervals r k = t k +i — Moreover, the low frequency 
noise is defined by the statistical properties of the signal at large-time-scale, i.e., by the 
fluctuations of the time difference A (k; q) at large q, determined by the slow dynamics 
of the average interpulse interval (q) = A(k;q)/q between the occurrence of pulses k 
and k + q. In such a case quite generally the dependence of the interpulse time r& on the 
occurrence number k may be described by the general Langevin equation with the drift 
coefficient d (r^) and a multiplicative noise h (r^) £ (k), 

d ^ = d{r k )+b{r k )ak). (11) 



Here we interpret k as continuous variable while the white Gaussian noise £ (k) satisfies the 
standard condition 

(t{k)t (*0> = *(*-*') 

with the brackets (...) denoting the averaging over the realizations of the process. Equation 
(II lj) we understand in Ito interpretation. 

Perturbative solution of Eq. (|11|) in the vicinity of r k yields 



k+j 

T k +j ^r k + d (r k ) j + b (r k ) J £ (I) dl, 

k 



k+q-1 



q k+j 



A (fc; g) = Yl Ti ~ J Tk +i dj - T k<l + d (r k ) y + b (r k ) J dj J £ (I) dl. 

l=k k 

After integration by parts we have 

2 k+Q 

A{k-q) = T k q + d{T k ) < ^ + b{T k ) J {k + q-l)£{l)dl, 



(12) 



(13) 



(14) 



(A(k;q)) = T k q + d(r k ) —. 



(15) 



Analogously, in the same approximation we can obtain and the variance a\(k;q) 
[A (k; q) 2 ) — (A (k; q)) 2 of the time difference A (k; q), 

a 2 A (k;q) = b 2 (r k ) —. 

A. Power spectral density 



(16) 



Substituting Eqs. (fT3j) and ([To]) into Eq. (fTU|) and replacing the averaging over k by the 
averaging over the distribution of the interpulse times r k we have the power spectrum 



S(f) = 4I 2 t / dr k P k (r k )Re / dqexpliuj 



r k q + d(r k ) 



Pk (T k ) Re 



e l v a 4 J er f c V /Z 



IX 



dr k 



(17) 



where x = irfr k /d(T k ). 

The replacement of the averaging over k and over realizations of the process by the 
averaging over the distribution of the interpulse times r k , P k (r k ), is possible when the process 
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is ergodic. Ergodicity is usually a common feature of the stationary process described by 
the general Langevin equation j^]. Therefore, we will consider the stationary processes 
of diffusion of the interpulse time T& described by Eq. (jlljl and restricted in the finite 
interval the motion. Such restrictions may be introduced as some additional conditions to 
the stochastic equation. The similar restrictions, however, may be fulfilled by introducing 
some additional terms into Eq. ([11)1 . corresponding to the diffusion in some "potential well", 
as in paper p. 

Approach (jl7j) is the improvement of the simplest model of the pure 1/f noise 
taking into account the second, drift, term d (r^) q 2 /2 in expression for A (k; q). Note, that 
for <i (Yfc) ^ from Eq. (|17|) we recover the known result 

According to Eqs. (JIJ, and (jl7j) the small interpulse times and the clustering of 
the pulses make the greatest contribution to 1/f 13 noise. The power-law spectral density is 
very often related with the power-law behavior of other characteristics of the signal, such as 
autocorrelation function, probabili ty d ensities and other statistics, and with the fractality 



of the signals, in general 
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35| . Therefore, we investigate the power-law 



dependences of the drift coefficient and of the distribution density on the time in some 
interval of the small interpulse times, i.e., 

d (r fc ) = jt^, Pk (r fc ) = Cr^ , r min < r k < r max (18) 

where the coefficient 7 represents the rate of the signal's nonlinear relaxation and C has to 
be defined from the normalization. 

The power-law distribution of the interpulse, interevent, interarrival, recurrence or wait- 
ing time is observable in different systems from physics, astronomy and seismology to the 



Internet, financial markets and neural spikes (see, e.g., |3j, ^ , [l5|, |36( and references herein). 

One of the most direct applications of the model described by Eq. (JT5j) . perhaps, is 
for the modeling of the computer network traffic ^| with the spreading of the packets of 
the requested files in the Internet traffic and exhibiting the power-law distribution of the 
inter-packet time. The modeling of these processes is under way. 

Because of the divergence of the power-law distribution and requirement of the stationar- 
ity of the process the stochastic diffusion may be realized over a certain range of the variable 
Tfc only. Therefore, we restrict the diffusion of 77 in the interval [T mm ,r max ] with the appro- 
priate boundary conditions. Then the steady state solution of the stationary Fokker-Planck 
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equation with a zero flow corresponding to Eq. (jllj) 



is Q 



C I „ ? d M 



For the particular power-law coefficients d{jk) and 6(r/c) (see, e.g., Eq. (J2SI)) we can obtain 
the power-law stationary distribution density (fTHj) . 

Then equations (fTTj) and (fTH|) yield the power spectra with different slopes /3, i.e., 

2 J 2 f fo\^ 
S (/) = ^/o V /k ( Xmin > Xmax ) ' ( 20 ) 



-Cmax 

4 (x m w,£ max ) = Re J e~^ x ~^ erfc (V-^c) x K rfx. (21) 



Here H o—X 2' "^rmn f / f'21 ■''max ,// 



fo = l(Cf)^, h = Z ^ =g , h = ~^s- (22) 



7T 7TTZ." TIT ■ 

max mm 



Note that fo is indefinite when a — > 0, however, / 2 " 4 is definite and converges to Cf in this 
limit. 

We note the special cases of the power spectral density (J2DJ). 
(i) h « / « / 2 , -K « < 1/2, 



r (1 + k) j //o 



5 (/) (2 - 5) cos [(*/2 + 1/4) tt] / V / ' ' 23 ' 



i.e., S (f) ~ 1// 1+2 -' 5 and S (/) ~ 1// for a = 0, in accordance with Eq. 
(ii) /</i,«> -1, 



s (fl-(rT^(f)"V57- (24) 



i.e., for very low frequencies 5 (/) ~ 1/V7- 
(iii) / » / 2 , « < 1/2, 

*</) = W(^)(!f F (25) 

i.e., for high frequencies 5 (/) ~ l// 2 - 

For very high frequencies / ^> t~2 x , however, we can not replace the summation in Eq. 
(ITUJ) by the integration. Then from Eqs. (jBJ) or (fTUj) one gets the shot noise S 1 (/) = 2al, 
Eq. 0. 
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Equations (f2T)J) and (|2*3|) - (j2l)|) reveal that the proposed model of the stochastic multiplica- 
tive point process may result in the power-law spectra over several decades of low frequencies 
with the slope f3 between 0.5 and 2. 

The simplest and well-known process generating the power-law probability distribution 
function for r k is a multiplicative stochastic process with b (r k ) = ar^ and 8 = 2/i — 1, 



written as 37] 



T fc+ i = r k + 7r fc 2M 1 + (JT%e k . (26) 

Here 7 represents the relaxation of the signal, while r k fluctuates due to the perturbation 
by normally distributed uncorrelated random variables e k with a zero expectation and unit 
variance and a is a standard deviation of the white noise. According to Eq. (JT§j) the steady 
state solution of the stationary Fokker-Planck equation with a zero flow, corresponding to 
Eq. (|26j) . gives the power- law probability density function for r k in the fc-space 

P*(T k )= l + _ a a = ^-2 M . (27) 

max min 



The power spectrum for the intermediate /, fi <C / <C f2, according to Eq. (j23j) is 

q( n (2 + «) (P-l)a 2 T Q9 -1/2) l_ 
UJ v^«(r^-r^)sin(7r/?/2)U //» 1 J 



where 



"= 1 + 3^' ^' <?<2 ' < 29) 



For /i = 1 we have a completely multiplicative point process when the stochastic change of 
the interpulse time is proportional to itself. Multiplicativity is an essential feature of the 
financial time series, economics, some natural and physical processes 

Another case of interest is with \x = 1/2, when the Langevin equation in the actual time 
takes the form 

%=T L + °Z (*) . (30) 
at r 

i.e., the Brownian motion of the interpulse time with the linear relaxation of the signal 
I ~ d/r. 

Figures [T] and EJ represent the spectral densities (jU with the different slopes (3 of the 
signals generated numerically according to Eqs. and (|26|) for the different parameters 
of the model. We see that the simple iterative equation (}2T)|) with the multiplicative noise 
produces the signals with the power spectral density of different slopes, depending on the 
parameters of the model. The agreement of the numerical results with the approximate 
theory is quite good. 
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FIG. 1: Power spectral density (@J) vs frequency of the signal generated by Eqs. © and (|26|) . 
Numerical simulations are averaged over 10 realizations of N = 10 6 pulse sequences with the 
parameters a = 1, [i = 1/2, a = 0.02 and different relaxations of the signal 7. We restrict the 
diffusion of the interpulse time in the interval r m i n = Tmax — 1 with the reflective boundary 

condition at r m \ n and transition to the white noise, t^+i = T max + ae^, for > r max . The straight 
lines represent the analytical results according to Eq. (|28p. 
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FIG. 2: The same as in Fig. ^but for fj, = 1, a = 0.05 and different parameters 7. 
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It should be noted that the low frequency noise is insensitive to the small additional 
fluctuations of the particular occurrence times t k . Therefore, we can interpret that Eqs. 
(f2U|) and (JHU|) describe the slow diffusive motion of the average interpulse time, superimposed 
by some additional randomness. 

On the other hand, the numerical investigations have shown that the proposed model is 
stable with respect to variation of dynamics of the interpulse time r k . The substitution of 
the reflecting boundaries for r k by an appropriate confining potential as in Ref. \^ does not 
change the result. 



B. Distribution density of the signal intensity 



The origin for appearance of 1/f fluctuations in the point process model described by Eqs. 
(|2*j)- (j3T))) is related with the slow, Brownian fluctuations of the interpulse time r k as a function 
of the pulse number k, when the average interpulse time Tfc(g) is a slowly fluctuating function 
of the arguments k and q. In such a case transition from the occurrence number k to the 
actual time t according to the relation dt = r k dk yields the probability distribution density 
Pt( T k) °f T k in the actual time t, 



Pt(r k ) = P k {T k )T k /'i 



(31) 



The signal averaged over the time interval r k according to Eq. (JHl) is I = ajr k . Therefore, 
the distribution density of the intensity of the point process (JBJ) averaged over the time 
interval r k is 

PC) - gft (?) ■ (32) 
If Pk{T~k) — const when r k — > (the condition for the exhibition for the pure 1/f noise in 
the point process model) the distribution density of the signal is 



P(I) 



(33) 



For the generalized multiplicative processes (jHJ), (JTTJ) , and |T8|) we have from Eqs. (J2 
and the distribution density of the signal intensity 

A- ] 



P{I) 



-A-1 



A-1 



T, 



A = 3 + a. 



(34) 



The power-law distribution of the signals is observab 



from earthquakes to the financial time series 
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e in a la rge variety of systems ranging 
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One of the simplest modelsgenerating the Brownian fluctuations of the interpulse time 
Tfc is an autoregressive model U Isl Q| with random increments and linear relaxation of the 
interpulse time, i.e., the model described by the iterative equation 

Tfc+i = T fc - 7(Y fe - f ) + ae k . (35) 

Here f is the average interpulse time, 7 is the rate of the linear relaxation, {eh} denotes the 
sequence of uncorrelated normally distributed random variables with zero expectation and 
a unit variance and a is the standard deviation of this white noise. The model (|3*|). (|1U|). 
and (J33j) then results in the power spectral density js| 



S(f) = P?f, a H = ^Ke~ K \ K = I^lI. (36) 



The distribution density of the intensity of the signal according to Eqs. (fT9")) and then 



is 

r2 f ^2 



Restricting the diffusion of the interpulse time 77 by the reflective boundary condition at 
T min > and for r min — > we have the truncated distribution density of the signal intensity 

P ' {,) ^/?ll + erf(A-)l eXP \ I 1 - 1 ) \W 1 > °' (38) 



/7F [1 + erf (AT)] 

In the asymptotic I I and / ^> 2K 2 I from Eq. (jHEJ) we have 




P r (I) ~ <^ ~ 1 (39) 
i.e., the power-law distribution density of the signal. Here 

a '» = TT^WY m 

The restriction of motion of by the reflective boundary condition at 77 = reduces 
the effective (average) value of Pk (0) = \ \P\. (rfc — > +0) + Pk (jk ~ > ~~ 0)] in Eq. (JTJ) and, 
consequently, the power spectral density approximately 2 times in comparison with the 
theoretical result ()36|) obtained without the restriction, because P& (r^ — > — 0) = for the 
restricted motion. More exactly, in such a case the power spectral density may be expressed 
by Eq. (J36j) with a r H instead of an, i.e., 

72 a H 



S r (f) = I 2 ^. (41) 
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C. Correlation function of the point process 



Correlation function C(s) of the point process (j3J) may be expressed as 

C ^ = ( ^ E -t k ~s)\=IaJ2 f (A (?)) S(A(q) - s)dA (q) = la ^ V q (s) 

\ k,q I 1 -oo 1 

(42) 

where the brackets (...) denote the averaging over the realizations of the process and over 
time (index k) as well. Such averaging coincides with the averaging over the distribution of 
the time difference A(q), ^ q (A (q)). 
From Eq. (}4~2*|) for the approximation 



k+q 

A(k; q) = t k+q -t k = J2 T i ~ r (<?)<?, ^ > ( 43 ) 

l=k+l 

we have the expression for the correlation function in the simplest approximation 



Tim ax 

C(s) ~IaJ2 Pk(r k )S(r k q - s)dr k = Ia5(s) + Ja^P fc (^ ' 

' mm 

Replacing the summation in Eq. (|44jl by the integration we have the approximate expression 
for the correlation function of the point processes (j3J) and (fTTj) or ([53]) 



(44) 



oo 

C(s) -Iajp k (~\ y, s > 0, C(-s) = C{s). 



(45) 



IV. SIGNAL AS A SUM OF UNCORRELATED COMPONENTS 

As it was already mentioned above, 1/f noise is often modeled as the sum of the Lorentzian 
spectra with the appropriate weights of a wide range distribution of the relaxation times 
r rel . It should be noted that the summation of the spectra is allowed only if the processes 
with different relaxation times are isolated one from another [fj, |3| . For the construction 
of the signal I(t) with 1/f noise spectrum from the stochastic equations with a wide range 
distribution of the relaxation times (and rates 7/ = l/r[ el ) one should express the signal as 
a sum of N uncorrelated components jj| 

N 

i( t ) = Y,m (46) 
1=1 

where every component l\ satisfies the stochastic differential equation 

ii = -7i{Ii-Ii) ( 4 7) 
13 



Here // is the average value of the signal component l\ , is the 5-correlated white noise, 
(£,i(t)£,i'(t')) — &i,v&{t — f), an d &i is the intensity (standard deviation) of the white noise. 
The distribution density P(Ii) of the component is Gaussian 

p ^-^7,^{-^ 1 '- 1 ^}- (48) 

The distribution density P(I) of the signal I(t), Eq. (jH3), expressed as a sum of uncor- 
rected Gaussian components, is Gaussian as well, 

p{1) = ^rM- [I if\' (49) 



with the average value / and the variance a 2 expressed as 



' = <* 2 = £ I:- (50) 



Therefore, the Bernamont-Surdin-McWhorter model based on the sum of signals with a 
wide range distribution of the relaxation times always results in the Gaussian distribution 
of the signal intensity. However, not all signals exhibiting 1/ / noise are Gaussian |2J]. Some 
of them are non-Gaussian, exhibiting power-law distribution or even fractal P, I30L 31, 32. 

HQS- 

Eqs. ({IE)) and ()47j) result in the expression for the correlation function of the signal (|4fij). 

C(s) = J2^ lS > S > Q - (51) 

The correlation function (JoTj) yields the power spectrum 

2rr 2 

Introducing the distribution of the relaxation rates, (7(7), we can replace the summation 
in Eqs. (I46J) and (J50j) - ([52j) by the integration and express the power spectrum of the signal 
(JIHJ) as 

7max £/max 

= / 2^ = 1 r a^p) 

J 7 + ^ ttj J 1 + y 

■Tmin 2/min 

Here 7 m i n and 7 max are minimal and maximal values of the relaxation rate, respectively. 
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A. Signals with the pure 1/ f power spectrum 



Eq. yields the pure 1/ / power spectrum only in the case when a 2 {uiy)g{uiy) = A 
const. In such the correlation function (JoTj) may be expressed as 



7max ' max 



CM = i J = f J (54) 

7min T r ^ 
min 

while the power spectrum (JHH)) yields 

#(/) = — 7 arctan _ arc t a n f^E^ ~ A 7min < w < 7max . (55) 

7T r L \ u / V w /J zf 



For the signal expressed not as a sum (J46j) but as an average of N uncorrelated compo- 
nents, 

1 N 

/ «( t ) = iv^ /i(t) ' (56) 

i=i 

all characteristics ([48 )1 - ()55 )) are similar, except that the average value I a of the averaged 
signal (jo^j) is N times smaller than that according to Eq. (|5U|) , while the expressions for the 
correlation function C(s), Eqs. (J5l"|) and (jUj), for the power spectrum S(f), Eqs. (}52|) . (|53|) . 
and (JHHJ) . and for the variance o^, Eq. fjBTljl . should be divided by N 2 , i.e., 

= °2 = ^£S- (57) 

- 2^/. ( 58 ) 

7max 

G»(«) = -^ I ^ 2 ( 7 )^( 7 )rf 7 . (59) 

Tin in 

When replacing the summation in Eqs. (j4"E|) . (|HH|l - (j5l^l and dSHJ) - (JHHJ) by the integration, 
the distribution density of the relaxation rates, g (7), should be normalized to the number 
of uncorrelated components N, 



<7( 7 )d 7 = N. (60) 



7 min 



We see the similarity of expressions (|45jl and ()59j1 for the correlation function of the point 
process model and that of the sum of signals with different relaxation rates, respectively. In 
general, however, different distributions Pk{Tk) of the interpulse time t*. when Pfc(0) 7^ 0, e.g., 
exponential, Gaussian and continuous distributions, with the slowly fluctuating interpulse 
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FIG. 3: Power spectra: (a) numerically calculated for the average signal ([5fi|l from N = 10 com- 
ponents ([47|) with J = 20, of( r yi)g('yi) = const and uniform distribution of lg7; of 7/ values in 
the interval 10~ 4 — 10°, i.e., with 3(7/) ~ jf , erf (7/) ~ 7;, and cri (71) = 0.1, open circles, in 
comparison with theoretical results (|58|) . straight line; (b) for the point process (jHJ), Q, and (|35|) 
with a = 1, f = 1, cr = 0.01, and 7 = 0.0001 averaged over 10 realizations of 10 5 pulse sequences, 
open circles, in comparison with the theoretical results according to Eq. (|41|). straight line, (c) 
and (d) numerically calculated distribution densities of the corresponding signals, open circles, in 
comparison with the theoretical results (|49|) . (|57|). and (|38|). solid lines, respectively. 

time may result in 1// noise. Therefore, the point process model is, in some sense, more 
general than the model based on the sum of the Lorentzian spectra. 

In Figure El the examples of the pure 1/ / power spectra for the average ()56|) of signals 
([47]) generated for different relaxation rates 7; and with the corresponding intensities of the 
white noise of and those of the autoregressive point process ©,(01), and ([Ho]) are presented 
together with the distribution densities of the corresponding signals. We see the similarity 
of the spectra but very different distributions of the intensity of the signals: the signal of 
the sum of the Lorentzians is Gaussian while that of the point process is approximately of 
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the power-law type, asymptotically P(T) ~ / 



B. Signals with the power spectral density of different slopes (3 



Using the sum of different Lorentzian signals we can generate not only a signal with the pure 
1/ / spectrum but the signal with any predefined slope (3 of power spectral density, as 
well. Indeed, let us investigate the case when 



(61) 



where A and rj are some parameters. Substitution of Eq. (joTjl into Eq. (J53J) yields the power 
spectral density 

s(/)= 4 / ^L dv 

TTj J i + y 

7min/w 



.4 



UJ 



1-7? 



7n 



1 



7n 



1 2 



7min 
Co> J 



,1, 



r/ + 1 



77+1 



7* 



r? + 1 



where s, a) is a Lerch's Phi transcendent. In the limit when 7 m ; n — > and 7„ 
we can approximate the power spectral density (jo^jl as 

o (n .. (g?Q!i 1 

UJ ~ 2cos(7rr 7 /2)/ 1 -'?' 

i.e., we have the generalization of the result (|55j) . 
For the average signal (jo^j) we have 

(2tt)M 1 



(62) 



■DC' 



(63) 



(64) 



2iV 2 COS (7T7//2) f 1 -^ 

In order to obtain an arbitrary (3 of the 1/ power spectral density we should choose in 
Eq. dSH v = 1 - P- 

The distribution density P a (I a ) of the average signal J a (t) is Gaussian 



Pa(/a) 



(J--fa) z 



2vro- a 



(65) 



with the variance a 2 expressed as 



2N 2 



E 

^ 2 (7k(7) d7 = ^(72 1 ax-7lJ 



7 



2iV 2 77 



(66) 
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The correlation function in such according to Eq. (J53|) is 

7max 

Ca(s) = 4f 2 \ e~^-i dl = _A_ [r (?7) 7minS) _ r ^ 7maxS)] 

Tmin 

where V (a, 2) is the incomplete gamma function. 
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FIG. 4: Power spectra: (a) numerically calculated for signal ()47|) . ()56j) and ()61j) from 10 components 
with I = 20, A = 100, 77 = —0.25, open circles, and r/ = 0.25, open squares, in comparison with 
theoretical results (|64|) . straight line; (b) for the point process ©, @ and l|26|) with the parameters 
a = 1, /i = 0.5, a = 0.02, and 7 = 0.0001, open squares, and 7 = 0.0003, open circles, averaged over 
10 realizations of 10 6 pulse sequences in comparison with the theoretical results (|28|) . straight lines, 
(c) and (d) numerically calculated distribution densities of the corresponding signals in comparison 
with the theoretical results (|65|) . (|66j) . and (|cS4l) . respectively, solid lines. 

Figure 0] demonstrates the possibility to generate stochastic signals exhibiting similar 1/ f 13 
power spectral densities with different slopes j3 by the summation of signals with different 
relaxation rates and according to the multiplicative point process model. The distribution 
densities of the corresponding signals are, however, completely different. 
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V. CONCLUSIONS 



The generalized multiplicative point processes (j2J), (JUJ), (jTHjl . and (J21j|) may generate time 
series exhibiting the power spectral density S(f) ~ l//' 3 with 0.5 < (3 < 2, Eqs. (fTTj) . 
(J2SJ) 7 an d (|25ji. i.e., with the slope observable in a large variety of systems. Such spectral 
density is caused by the stochastic diffusion of the interpulse time, resulting in the power-law 
distribution. The power-law distribution of the interpulse, interevent, interarrival, recurrence 
or waiting times is observed in different systems from physics, astronomy and seismology to 
the Internet, financial markets, neural spikes, and human cognition. 

Furthermore, the power-law distribution of the interpulse time results in the power-law 
distribution of the stochastic signal, P(I) ~ I~ x with 2 < A < 4, i.e., the phenomenon 
observable in a large variety of processes, from earthquakes to the financial time series, as 
well. The proposed model relates and connects the power-law autocorrelation and spectral 
density with the power-law distribution of the signal intensity into the consistent theoretical 
approach. The generated time series of the model are fractal since they exhibit jointly the 
power-law probability distribution and the power-law autocorrelation of the signal. 

In addition, we have analyzed the relation of the point process model with the Bernamont- 
Surdin-McWhorter model of 1/f noise, representing the signal as a sum of the appropriate 
signals with the different rates of the linear relaxation. From the performed analysis we can 
conclude that the multiplicative point process model of 1/f noise when the signal consisting 
of pulses with a stochastic motion of the interpulse time is more general and complementary 
to the model based on the sum of signals with a wide-range distribution of the relaxation 
times. In contrast to the Gaussian distribution of the intensity of sum of the uncorrelated 
components, the point process model generating 1/f noise exhibits the power-law distribu- 
tion of the intensity of the signal. Moreover, it is free from the requirement of a wide-range 
distribution of the relaxation times. Obviously, the multiplicative point process model of 
1/f 13 noise may be used for modeling and analysis of stochastic processes in different systems 
exhibiting the pulsing signals. 
Acknowledgments 

We acknowledge support by the Lithuanian State and Studies Foundation. 

[1] J. B. Johnson, Phys. Rev. 26, 71 (1925); W. Schottky, Phys. Rev. 28, 74 (1926); D. W. Allan, 
Proc. IEEE 54, 221 (1966); C. W. J. Granger, Econometrica 34, 150 (1966); R. F. Voss and J. 

19 



Clarke, Nature (London) 258, 317 (1975); W. H. Press, Comments Astrophys. 7, 103 (1978); 
M. Gardner, Sci Am. 238, 16 (1978); A. Lawrence, M. G. Watson, K. A. Pounds, and M. 
Elvis, Nature (London) 325, 694 (1987). 

[2] F. N. Hooge, T. G. M. Kleinpenning and L. K. J. Vadamme, Rep. Prog. Phys. 44, 479 (1981); 
P. Dutta and P. M. Horn, Rev. Mod. Phys. 53, 497 (1981); S. M. Kogan, Usp. Fiz. Nauk 
145, 285 (1985) [Sov. Phys. Usp. 28, 170 (1985)]; M. B. Weissman, Rev. Mod. Phys. 60, 537 
(1988); C. M. Van Vliet, Solid State Electron. 34, 1 (1991); G. P. Zhigal'skii, Usp. Fiz. Nauk 
167, 623 (1997) [Phys.-Usp. 40, 599 (1997)]; E. Milotti, arXiv: |physics/0204"033| (2002); H. 
Wong, Microel. Reliab. 43, 585 (2003). 

[3] S. Thurner, S. B. Lowen, M. C. Feurstein, C. Heneghan, H. G. Feichtinger, and M. C. Teich, 
Fractals 5, 565 (1997). 

[4] D. L. Gilden, T. Thornton, and M. W. Mallon, Science 267, 1837 (1995). 

[5] B. Kaulakys and J. Ruseckas, Phys. Rev. E 70, 020101(R) (2004); |cond-mat/0408507l 

[6] J. Bernamont, Ann. Physik 7, 71 (1937); M. Surdin, J. Phys. Radium (Serie 7) 10, 188 
(1939); F. K. du Pre, Phys. Rev. 78, 615 (1950); A. Van der Ziel, Physica 16, 359 (1950); A. 
L. McWhorter, in Semiconductor Surface Physics, ed. R. H. Kingston, (University of Penn- 
sylvania, Philadelphia, 1957), p. 207; V. Palenskis, Lit. Fiz. Sb. 30, 107 (1990) [Lith. Phys. J. 
30(2), 1 (1990)]; F. N. Hooge, in Proc. lJ^th Int. Conf. on Noise in Physical Systems and 1/f 
Fluctuations, Leuven, Belgium, 14-18 July 1997, edited by C. Claeys and E. Simoen (World 
Scientific, Singapore, 1997), p. 3. 

[7] B. Kaulakys, arXiv: adap-org/9806004 (1998). 



[8] B. Kaulakys and T. Meskauskas, Phys. Rev. E 58, 7013 (1998); adap-org/9812003; B. 



Kaulakys, Phys. Lett. A 257, 37 (1999); adap-org/9907008 



[9] B. Kaulakys, Lith. J. Phys. 40, 281 (2000). 
[10] B. Kaulakys and T. Meskauskas, Nonlin. Anal.: Mod. Contr. (Vilnius) 4, 87 (1999); B. 

Kaulakys and T. Meskauskas, Microel. Reliab. 40, 1781 (2000); |cond-mat/0 303603 B. 

Kaulakys, Microel. Reliab. 40, 1787 (2000); |cond-mat/0305067l 
[11] K. L. Schick and A. A. Verveen, Nature (London) 251, 599 (1974); F. Gruneis, Physica A 

123, 149 (1984); F. Gruneis and H.-J. Baiter, Physica A 136, 432 (1986); F. Gruneis and T. 

Musha, Jpn. J. Appl. Phys. 25, 1504 (1986); T. Musha, Jpn. J. Appl. Phys. 26, 2022 (1987); 

F. Gruneis, M. Nakao, and M. Yamamoto, Biolog. Cybern. 62, 407 (1990); F. Gruneis, M. 

Nakao, Y. Mizutani, M. Yamamoto, M. Meesmann, and T. Musha, Biolog. Cybern. 68, 193 

(1993); F. Gruneis, Fluct. Noise Lett. 1, R119 (2001); F. Gruneis, Fluct. Noise Lett. 4, L413 

(2004). 



20 



[12] H. R. Bittner, P. Tosi, C. Braun, M. Meesmann, and K. D. Kniffki, Geolog. Rundschau 85, 
110 (1996); R. Console and M. Murru, J. Geophys. Res. 106, 8699 (2001); L. Telesca, V. 
Cuomo, and V. Lapenna, Fluct. Noise Lett. 2, L357 (2002); J. M. Halley and P. Inchausti, 
Fluct. Noise Lett. 4, Rl (2004); M. Lindman, K. Jonsdottir, R. Roberts, B. Lund, and R. 
Bodvarsson, Phys. Rev. Lett. 94, 108501 (2005). 

[13] T. Musha and H. Higuchi, Jap. J. Appl. Phys. 15, 1271 (1976); X. Zhang, G. Hu, Phys. Rev. 
E 52, 4664 (1995); M.Y. Choi and H.Y. Lee , Phys. Rev. E 52, 5979 (1995). 

[14] I. Csabai, J. Phys. A 27, L417 (1994); W. Willinger, M. S. Taqqu, R. Sherman, and D. V. 
Wilson, IEEE/ACM Trans. Networking, 5, 71 (1997); I. Kaj and I. Marsh, In: Quality of 
Service in Multiservice IP Networks, ed. M. Ajmone Marsan, G. Corazza. M. Listanti, A. 
Roveri (Springer- Verlag, Heidelberg, 2003), Lecture Notes in Computer Science 2601, p. 35; 
A. J. Field, U. Harder, and P. G. Harrison, IEE Proceedings-Communications 151, 355 (2004); 
S. Cohen, M. S. Taqqu, Methodol. Comput. Appl. Probabil. 6 363 (2004); L. A. N. Amaral, 
A. Diaz-Guilera, A. A. Moreira, A. L. Goldberger, and L. A. Lipsitz, Proc. Natl. Acad. Sci. 
USA 101, 15551 (2004). 

[15] J. Davidsen and H. G. Schuster, Phys. Rev. E 62, 6111 (2000); J. Davidsen and M. Paczuski, 
Phys. Rev. E 66, 050101(R) (2002); Y. Kamitani, H. Suyari, and I. Matsuba, Electron. Com- 
munic. Jap. Pt.III 87, 72 (2004). 

[16] J. I. Yano, K. Fraedrich, and R. Blender, J. Climate 14, 3608 (2001); W Muller, R. Blender, 
and K. Fraedrich, Nonlin. Process. Geophys. 9, 37 (2002); W. W. Tung, M. W. Moncrieff, and 
J. B. Gao, J. Climate 17, 2736 (2004); J. I. Yano, R. Blender, C. D. Zhang, and K. Fraedrich, 
Quart. J. Royal Meteor. Soc. 130, 1697 (2004); K. Fraedrich, U. Luksch, and R. Blender, 
Phys. Rev. E 70, 037301 (2004). 

[17] C. Guidorzi, F. Frontera, E. Montanari, M. Feroci, L. Amati, E. Costa, and M. Orlandini, 
Astron. Astrophys. 416, 297 (2004). 

[18] M. Agu, H. Akabane, and T. Saito, J. Phys. Soc. Jap. 70, 2798 (2001). 

[19] M. Celasco and R. Eggenhoffner, Eur. Phys. J. 23, 415 (2001). 

[20] J. Davidsen and H. G. Schuster, Phys. Rev. E 65, 026120 (2002). 

[21] V. Gontis, Lith. J. Phys. 41, 551 (2001); arXiv: |cond-mat/0201514| (2002); V. Gon- 
tis and B. Kaulakys, Proc. 17th Intern. Conf. on Noise and Fluctuations, edited by J. 
Sikula (Prague, August 18-22, 2003) p. 622; V. Gontis, B. Kaulakys, M. Alaburda, and 
J. Ruseckas, Sol. St. Phenom. 97-98, 65 (2004); V. Gontis and B. Kaulakys, Physica A 
343, 505 (2004); |cond-mat/0303089j V. Gontis and B. Kaulakys, Physica A 344, 128 (2004); 



cond-mat/0412723 



21 



[22] A. Coza and V. V. Morariu, Physica A 320, 449 (2003); E. -J. Wagenmakers, S. Farrell, and 

R. Ratcliff, Psychonomic Bulletin & Review 11, 579 (2004). 
[23] J. B. Gao, Medic, k Biolog. Engineer. & Comput. 42, 345 (2004). 
[24] D. Delignieres, L. Lemoine, and K. Torre, Human Movement Science 23, 87 (2004). 
[25] M. B. Priestley, Spectral Analysis and Time Series (Academic Press, San Diego, 1981). 
[26] J. Timmer and M. Konig, Astron. Astrophys. 300, 707 (1995). 

[27] S. L. Marpl, Jr., Digital Spectral Analysis with Applications (Prentice-Hall, New Jersey, 1987). 
[28] P. Stoica and R. L. Moses, Introduction to Spectral Analysis (Prentice-Hall, New Jersey, 1997). 
[29] C.W. Gardiner, Handbook of Stochastic Methods for Physics, Chemistry and Natural Sciences 

( Springer- Verlag, Berlin, 1985). 
[30] B. B. Mandelbrot, J. Business 36, 394 (1963); I. McHardy and B. Czerny, Nature (London) 

325, 696 (1987); G. W. Wornell, Proc. IEEE 81, 1428 (1993); H. Scher, G. Margolin, R. 

Metzler, J. Klafter, and B. Berkowitz, Geophys. Res. Lett. 29, 1061 (2002); P. Bak, K. 

Christensen, L. Danon, and T. Scanlon, Phys. Rev. Lett. 88, 178501 (2002); L. Giraitis, P. 

Kokoszka, R. Leipus, and G. Teyssiere, Acta Applicandae Mathematicae 78, 285 (2003). 
[31] B. B. Mandelbrot, Fractals and Scaling in Finance (Springer, New York 1997). 
[32] B. B. Mandelbrot, Multifractals and 1/f Noise: Wild Self- Affinity in Physics (Springer, New 

York 1999). 

[33] R. N. Mantegna, H. E. Stanley, An Introduction to Econophysics: Correlations and Complexity 

in Finance (Cambridge University Press, Cambridge, 2000). 
[34] S. T. Rachev (Editor), Handbook of Heavy Tailed Distributions in Finance (Elsevier B.V., 

2003). 

[35] X. Gabaix, P. Gopikrishnan, V. Plerou, and E. Stanley, Nature (London) 423, 267 (2003). 
[36] S. B. Lowen and M. C. Teich, IEEE Trans. Inform. Theory 39, 1669 (1993); M. Usher, M. 

Stemmler, and Z. Olami, Phys. Rev. Lett. 74, 326 (1995); M. Takayasu, H. Takayasu, and 
T. Sato, Physica A 233, 824 (1996); M. S. Wheatland, Astrophys. J. 536, L109 (2000); J. 
W. Kirchner, X. Feng, and C. Neal, Nature (London) 403, 524 (2000); E. Scalas, Physica A 
287, 468 (2000); F. Lepreti, V. Carbone, and P. Veltri, Astrophys. J. 555, L133 (2001); P. 
Richmond, Eur. Phys. J. 20, 523 (2001); L. Sabatelli, S. Keating, J. Dudley, and P. Richmond, 
Eur. Phys. J. 27, 273 (2002); A. Erramilli, M. Roughan, D. Veitch, W. Willinger, Proc. IEEE 
90, 800 (2002); S. Abe and N. Suzuki, Eurohys. Lett. 61, 852 (2003); P. L. Conti, A. Lijoi, 
and F. Ruggeri, Appl. Stochastic Models Bus. Ind. 20, 305 (2004); F. Selcuk, Physica A 333, 
306 (2004); T. Kaizoji and M. Kaizoji, Physica A 336, 563 (2004); A. Corral, Phys. Rev. Lett. 
92, 108501 (2004); A.-H. Sato, Phys. Rev. E 69, 047101 (2004); M. E. J. Newman, arXiv: 



22 



|cond-mat/0412004| S. M. D. Queiros, C. Anteneodo, C. Tsallis, arXiv: physics/0503024, 
[37] V. Gontis, Nonlin. Anal.: Mod. Contr. (Vilnius) 7, 43 (2002); arXiv: |cond-mat/0211317| 
(2002). 

[38] A.-H. Sato, H. Takayasu, and Y. Sawada, Phys. Rev. E 61, 1081 (2000); Z. F. Huang and S. 

Solomon, Physica A 306, 412 (2002). 
[39] P. Gopikrishnan, M. Meyer, L. A. N. Amaral, and H. E. Stanley, Eur. Phys. J. B 3, 139 

(1998); W. J. Reed and B. D. Hughes, Phys. Rev. E 66, 067103 (2002); N. Scafetta and B. 

J. West, Phys. Rev. Lett. 92, 138501 (2004); C. B. Yang, J. Phys. A: Math. Gen. 37, L523 

(2004); L. A. N. Amaral, D. J. B. Soares, L. R. da Silva, L. S. Lucena, M. Saito, H. Kumano, 

N. Aoyagi, and Y. Yamamoto, Europhys. Lett. 66, 448 (2004); A. Fujihara, T. Ohtsuki, and 

H. Yamamoto, Phys. Rev. E 70, 031106 (2004). 
[40] F. N. Hooge, IEEE Trans. Electron. Dev. 41, 1926 (1994); F. N. Hooge, Physica B 239, 223 

(1997). 



23 



